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ABSTRACT 


Plecoglossus altivelis (ayu) is an amphidromous fish widely distributed in Northeastern Asia from the East 
China Sea to the northern Japanese coastal waters, encompassing the Korean Peninsula within its range. The 
shore lines of northeastern region in Asia have severely fluctuated following glaciations in the Quaternary. In 
the present study, we investigate the population genetic structure and historical demographic change of P. 
altivelis at a population level in East Asia. Analysis of molecular variance (AMOVA) based on 244 mitochon- 
drial control region DNA sequences clearly showed that as the sampling scope extended to a larger geographic 
area, genetic differentiation began to become significant, particularly among Northeastern populations. A 
series of hierarchical AMOVA could detect the genetic relationship of three closely located islands between 
Korea and Japan that might have been tightly connected by the regional Tsushima current. Neutrality and 
mismatch distribution analyses revealed a strong signature of a recent population expansion of P. altivelis in 
East Asia, estimated at 126 to 391 thousand years ago during the late Pleistocene. Therefore it suggests that 
the present population of P. altivelis traces back to its approximate demographic change long before the last 
glacial maximum. This contrasts our a priori expectation that the most recent glacial event might have the 
most crucial effect on the present day demography of marine organisms through bottleneck and subsequent 
increase of effective population size in this region. 
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lation demographic expansion 


INTRODUCTION 


Fluctuations of climate and sea level during the Quaternary 
ice ages have profoundly influenced the present population 
genetic structures of marine and terrestrial species (Hewitt, 
1996, 2000; Palumbi, 1997). Similarly associated with pale- 
oceanographic changes, a rapid population increase has been 
repeatedly reported in pelagic fishes such as tuna, angelfish, 
sardine and anchovy during the Pleistocene (Bowen and 
Grant, 1997; Lecomte et al., 2004; Alvarado Bremer et al., 
2005). In the Indo/West Pacific region, past climate changes 
are shown to have affected the evolutionary processes of 
speciation and population demographic changes for several 


teleost marine fishes (Palumbi, 1997; Williams and Benzie, 
1997, 1998; Duke et al., 1998; Briggs, 1999; Sugama et al., 
2002; Williams et al., 2002; Lourie and Vincent, 2004). 
Nevertheless, the effects of past climate and oceanographic 
changes on marine species in Northeast Asia remain relative- 
ly unknown, and thus interests in this question have recently 
began to grow. 

The Western Pacific marginal seas in contiguity with the 
Korean Peninsula, the East China Sea and the East/Japan 
Sea, are interesting from a paleoceanography perspective 
because they experienced extreme topographical changes 
during the late Quaternary period (Fig. 1B) (Wang, 1999). 
When the global sea level reached approximately 120-130 
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m lower than the present level during the last glacial maxi- 
mum (LGM; 19-26 thousand years ago) (Clark et al., 2009), 
the East/Japan Sea was closed by the rise of continental shelf 
that connected Eastern China with the southern part of the 
Japanese archipelago and the Yellow Sea had almost disap- 
peared due to shallowness (Xu and Oda, 1999; Myers et al., 
2000; Kong et al., 2006; Xiang et al., 2008). The present 
topography of this region was formed since sea level had 
risen during the post-glacial period and the Tsushima Current 
had flowed into the East/Japan Sea (Kitamura et al., 1999, 
2001; Ujiié and Ujiié, 1999; Park et al., 2000; Li et al., 2007). 
Accordingly, one may expect marine organisms in this region 
to have a general phylogeographic signature. For example, 
historical demographic expansion to shallow continental 
margins via habitat expansion due to a rising sea level and 
population genetic structures reflecting historical land bridges 
are the most expected trends. In fact, evidence of rapid popu- 
lation demographic expansion has been gathered from a few 
taxonomically distant marine species in East Asia, such as 
the Japanese anchovy, Engraulis japonicas (Liu et al., 2006), 
crimson snapper, Lutijanus erythropterus (Zhang et al., 2006), 
and sea slater, Ligia sp. (Jung et al., 2008). However, sam- 
ples were limited most of those studies, as they considered 
only a part of the ranges of the species within the region. 
Therefore, attempts to investigate an entire evolutionary 
history of one species which had been influenced from the 
paleoceanographic changes in the region have been hindered. 
Combining and comparatively analyzing molecular genetic 
data from multiple countries, which were reported by dif- 
ferent researchers, suggested to be helpful for solve these 
problems. 

In the present study, we demonstrate this approach through 
phylogeographic study of the ayu, Plecoglossus altivelis, in 
East Asia. The ayu, Plecoglossus altivelis (Order Osmeri- 
formes, family Plecoglossidae), is an osmerid fish that mig- 
rates between freshwater and the sea during juvenile period. 
This species is widely distributed in East Asia, especially in 
Korea, Japan, Taiwan and the entire Chinese coastal Provin- 
ces except Jiangsu. It is an interesting species to investigate 
the effects of past climate changes on the population demo- 
graphy, change of distribution range and thus molecular evo- 
lution of marine inhabitants in north East Asia during the late 
Pleistocene period because its distributional range almost 
overlaps with the seas surrounding the Korean Peninsula, 
where severe climate change would have left salient signa- 
tures on the evolution of marine animals. Previous studies 
on this species have focused on the ecological aspects of its 
amphidromous life history and on stock enhancement due to 
its high commercial values in Korea and Japan (Chyung, 
1977; Suda and Gomyo, 1995; Iguchi, 1996; Takahashi et al., 
1998; Azuma et al., 2003; Kim and Park, 2005; Yamamoto 
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et al., 2007). A few population genetic studies based on mol- 
ecular DNA sequence data have been carried out around the 
Japanese archipelagos, with an emphasis on population dif- 
ferentiation between the mainland and island populations 
(Iguchi, 1996; Iguchi et al., 2006). Therefore, to understand 
the phylogeography and effects of past climate changes on 
this species at the regional level, more extensive samples from 
other regions which considered being their habitats need to be 
analyzed. 

In this study, we used previously published data from the 
Japanese archipelagos combining our new mitochondrial 
control region DNA data. By doing this we intend to address 
the question of how and to what extent paleoceanographic 
changes during the Quaternary period would have impacted 
the population genetic structure and historical population 
demography of marine species in Northeast Asia. 


MATERIALS AND METHODS 


Sample collection 

Sixty-eight ayu individuals were caught from seven localities 
(Table 1) in Korea from March to July 2007, when the fishes 
were returning to rivers. The Korean mainland samples were 
acquired from six localities (Myungpa, Namdae, Osip, Maeup, 
Jusu and Wangpi streams) on the eastern coast and a single 
representative locality from the southern coast, the Sumjin 
River (SUM). These localities were specifically chosen to 
assess the fine scale genetic structure of ayu populations in 
Korea (Fig. 1B). 

To analyze population genetic structure of ayu fishes at 
an extended geographic scale of East Asia, 176 additional 
control regions of mtDNA sequences of Korean, Japanese 
and Chinese P. altivelis were acquired from previous stud- 
ies (Iguchi et al., 1999; Iguchi and Nishida, 2000) and Gen- 
Bank. Ten sequences from the Jejudo Island (JEJ), which is 
located in the Southern Sea of Korea, and 161 sequences from 
13 Japanese ayu populations were generously provided by 
Iguchi and Nishida upon our request. In addition, we used 
five sequences from a Chinese population, although we were 
unfortunately not able to obtain location information for the 
five samples (accession no. EU124679 -83). 


DNA extraction, amplification and sequencing 

Genomic DNA was isolated from pectoral fins using Labo- 
Pass Tissue Mini kits (Cosmo Genetech, Seoul, Korea). A 
partial fragment of about 430 base pairs (bp) from the mito- 
chondrial control region was amplified by polymerase chain 
reaction (PCR) using the primers L15923 (5- TTA AAG САТ 
CGG TCT TGT AA-3’) and H16498 (5’-CCT GAA GTA 
GGA ACC AGA TG-3’) (Shields and Kocher, 1991). The 
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Fig. 1. Haplotype network of mtDNA control region sequences and sample collection sites of Plecoglossus altivelis. А, The large and 
small circles, both called nodes, represent each unique haplotype observed and those unseen in the DNA samples, respectively. The 
black and white color with a different hatching pattern for the nodes indicates eight haplogroups, which can be grouped by their 
closeness in genetic distance to each other. Although genetic closeness is the primary basis for the grouping, we arbitrarily decided the 
membership of any ambiguous nodes into certain groups as shown in the figure for the sake of convenience in labeling. Two of the 
haplotypes in the box could not be connected to any other haplotype in 9596 parsimonious criterion, and they were from SAD and 
SAM respectively; B, Map of sampling localities and proportion of the mtDNA haplotypes defined in A. The darkly shaded sections on 
the map represent the current terrestrial area, and the lightly shaded area and its boundary lines represent ancient land masses and 
coastal lines during the last glacial maximum period (Ujiié and Ujiié, 1999). Each locality is represented by a pie chart and abbrevia- 
tion that correspond to Table 1. The pie charts display the proportion of clade defined in the haplotype network at the different 
sample site. 
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Table 1. Samples of ayu (Plecoglossus altivelis) from different localities in Korea, Japan, and China 


Sample site Abbreviation Type Ni Nu N/N; h+s.d. m+s.d. 
Korea (8 populations) 
Myungpa River MYUNG M 11 11 9 1.00 1.00--0.0388 0.0260--0.0147 
Namdae River NAM M 10 10 5 1.00 1.00--0.0447 0.0280--0.0160 
Jusu River JU M 10 10 7 1.00 1.00--0.0447 0.0265--0.0154 
Osip River OS M 10 8 4 0.80 0.96+0.0594 0.0241+0.0139 
Maeup River MA M 13 13 11 1.00 1.00--0.0302 0.0225--0.0126 
Мапарі Кіуег WANG M 4 4 3 1.00 1.00 0.1768 0.0263 0.0184 
Sumjin River SUM M 10 10 5 1.00 1.00+0.0447 0.0299+0.0169 
Jejudo Island JEJ I 10 5 3 0.50 0.76+0.1295 0.0078 +0.0051 
Ѕит 0.99+0.0050 0.0255+0.0330 
Japan (13 populations) 
Fukue Island FKE I 10 10 6 1.00 1.00--0.0447 0.0180--0.0106 
Tsushima Island TUS I 10 8 4 0.80 0.96+0.0594 0.0159+0.0095 
Chikugo River KYS M 15 15 11 1.00 1.00--0.0243 0.0296--0.0161 
Oki Island OKI I 10 3 2 0.30 0.38+0.1813 0.0024+0.0022 
Biwa Lake BWL L 15 5 4 0:33 0.56+0.1434 0.0036+0.0027 
Sado Island SAD I 10 10 9 1.00 1.00--0.0447 0.0282--0.0160 
Shinano River HNS M 15 15 12 1.00 1.00 0.0243 0.0299 0.0163 
Iwaki River IWK M 15 13 10 0.87 0.98+0.0308 0.0253+0.0139 
Same River SAM M 15 15 14 1.00 1.00--0.0243 0.0334--0.0180 
Kanzawa River KZW M 15 15 13 1.00 1.00 0.0243 0.0254 0.0140 
Amori River AMR M 15 14 14 0.93 0.99 0.0281 0.0247 0.0136 
Yaku Island YAK I 6 4 4 0.67 0.802 0.1721 0.01882 0.0120 
Amamioshima Island AMM I 10 4 4 0.40 0.80+0.0888 0.0092+0.0059 
Sum 0.99+0.0037 0.0280+0.0143 
China (1 population)? 
China CHI 5 4 4 0.80 0.90+0.1610 0.0036+0.0032 
Total 244 206 176 0.72 0.99+0.0021 0.0284+0.0007 


Type, sample site (M, mainland; I, island; and L, lake); М, number of individuals; Na, number of haplotypes; Nu, number of unique haplotypes; h, haplo- 


type diversity; and z, nucleotide diversity per site. 


?The information on the sampling sites of the Chinese samples is not known, but only five sequences were retrievable from GenBank (EU124679-83). 


20 uL PCR contained 2.0 uL 10x PCR buffer, 0.05 umol 
Mg?*, 0.001 umol for each dNTP, 8 pmol of each primer, 1 
ug bovine serum albumin, 0.3 unit taq polymerase and 50 ng 
genomic DNA of each sample. The temperature profile of 
the PCR was set at 94°C for 2 min, followed by 35 cycles of 
94°C for 1 min, 50°C for 45s and 72°C for 2 min 10s, with 
a terminal elongation at 72°C for 7 min. PCR products were 
purified with a LaboPass PCR purification kit (Cosmo Gene- 
tech). Sequencing was bidirectionally performed with an 
ABI PRISM 3100 automated sequencer (Applied Biosystems, 
Foster City, CA, USA). 


DNA analysis 

The 68 mtDNA control region sequences were checked and 
confirmed visually with the software AlignIR (Licor Co., Lin- 
coln, NE, USA). A total of 244 control region sequences, 
including 176 additional sequences from East Asia popula- 
tions, were aligned using clustalW (Thompson et al., 1994), 
which is implemented in the software MEGA (Tamura et al., 
2007). Finally, the DNA fragments were confined to a length 
of 329 bp for further comparisons and analyses (GenBank 
accession no. GQ409863-932). 
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The 68 DNA control region sequences were used in the 
analysis of the population genetics of Korean mainland P. 
altivelis populations as well as in the comparative analysis 
with the other Asian populations. Analyses of molecular vari- 
ance (AMOV A) were conducted by the population genetics 
computer program Arlequin version 3.11 (Excoffier et al., 
2005). To evaluate polymorphism of the control region DNA 
sequences, haplotype diversity (Л) (Nei, 1987) and nucleo- 
tide diversity (z) (Nei and Li, 1979) were calculated with 
Arlequin for all East Asian populations. 

A series of hierarchical AMOVA and pairwise Ф statistics 
were also conducted with the same program to draw a natur- 
al and regional boundary of P. altivelis population genetic 
structure in East Asia. Grouping for AMOVA was based on 
G1 value (Table 2). For some analyses, data from the Kor- 
ean mainland ayu populations were pooled into one Korean 
P. altivelis population because no genetic differentiation was 
found among them. Of the Japanese populations, the land- 
locked Biwa Lake population was excluded because of its 
high genetic differentiation with other populations in the pre- 
vious study (Iguchi et al., 1999). 

A parsimonious network of the haplotypes was construct- 
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ed with TCS (Clement et al., 2000). We visually divided the 
haplotypes in the minimum-spanning network into eight hap- 
logroups (Fig. 1A) according to an ad hoc criterion, approxi- 
mate genetic closeness among them, and then presented their 
frequency as pie diagrams for solely temporary comparisons 
between localities (Fig. 1B). 

To investigate population demographic changes in the re- 
cent past, the entire East Asian data set and each of the genet- 
ically distinct subsets were examined for neutrality of mark- 
ers with Tajima's D test (Tajima, 1989, 1993) and Fu's Fs 
test (Fu, 1997) using Arlequin (Excoffier et al., 2005). Addi- 
tionally, the mismatch distribution analysis (Schneider et al., 
1999) was conducted to test the model of sudden population 
expansion (Rogers, 1995) for the genetically distinct groups 
from the result of hierarchical AMOVA and whole East Asia 
data respectively. The genetic groups that consist of only one 
population were simply excluded because they are inappro- 
priate to represent population expansion by the effects of 
paleoceanographic changes in East Asia. 


RESULTS 


We obtained 244 control region sequences (329 bp) of P. alti- 
velis from 21 localities in East Asia (Korea, Japan, and China), 
including 176 unique haplotypes. From these samples, the 
total nucleotide diversity was 0.0284 + 0.0007, and the value 
of nucleotide diversity within each country was also low 
(Korea, 0.0255; Japan, 0.0280; and China, 0.0036) (Table 1). 
Within the Korean mainland, 68 mtDNA control region 
DNA sequences from seven localities were first examined 
for population genetic structure. AMOVA could not find any 
evidence of genetic differentiation, at either regional (i.e., 
eastern vs. southern coast) or population level, indicating 
high genetic homogeneity of P. altivelis populations in Korea. 
When the populations were examined according to two hypo- 
thetical groups, that is, Eastern and Southern (according to 
the geography of the Korean Peninsula) (Fig. 1B), the Фст 
value was estimated at —0.0277 (p=1.00). The great majo- 
rity of variation was found within populations (103.5696) 
rather than among groups (—2.77%) or among populations 
within groups (—0.79%). By contrast, when the JEJ data 
(Iguchi et al., 1999) was combined with the above Korean 
mainland data, highly significant and large pairwise Ф;т val- 
ues appeared for all the comparisons between JEJ and the 
other mainland populations (Table 2), indicating allopatric 
isolation of the JEJ insular population from the mainland. 
The pairwise Фет values among populations of ауа from 
the entire sampling region are shown in Table 2. We could 
find the following general patterns from these values: 1) 
almost all the pairwise Фот values between mainland popu- 


284 


lations within both Korea and Japan were low and not signi- 
ficant, being below the overall Фот value (0.157); ii) most 
pairwise Фот values between insular populations and main- 
land populations were significant except for those of Sado 
Island (SAD); iii) comparisons between insular populations 
resulted in high and significant Фет values; and iv) the iso- 
lated Biwa lake population also showed high and significant 
Фет values similar to the other insular populations (Table 2). 
However, the significances of Фет values from the most pair- 
wise comparisons disappeared after Bonferroni correction 
(p < 0.0001), except for the comparisons including Jejudo 
island and China populations. 

According to the hierarchical AMOVA based on the pair- 
wise Фет value, we could obtain a series of Фст values cor- 
responding to each grouping method. The highest Фст value, 
which indicates the highest proportion of genetic variance 
able to be attributed to groups given the present data, was 
observed as 0.210 when the populations were grouped into 
following seven groups: Korean mainland populations (refer- 
red as “Korea” hereafter); JEJ, Tsushima Island (TUS), Fukue 
Island (ЕКЕ) (Tushima"); Chikugo River (KYS), Shinano 
River (HNS), Iwaki River (IWK), Amori River (AMR), Kan- 
zawa River (KZW), Same River (SAM), SAD ("Japan"); 
Amamioshima Island (AMM); Yaku Island (YAK); Oki 
Island (OKI); and China. The results indicated that the Фст 
values tended to be higher when the three island populations 
(JEJ, TUS, FKE) located between the Korean Peninsula and 
Japanese archipelago were grouped together (Appendix 1). 
Of four distant Japanese islands, AMM, YAK, and OKI form- 
ed separate groups, consistent with high Фет values in Table 
2 and only SAD was included in the group of Japanese main- 
land populations. AMOVA of this grouping resulted in a 
significant hierarchical genetic structure in the ayu popula- 
tions (Ф;т=0.246, р< 0.000; Фст=0.210, p=0.000). The 
great majority of variation was contained within populations 
(75.40%) rather than among groups (21.02%) or among 
populations within groups (3.58%) (Table 3). 

The assumption of neutrality of the control region marker 
was evaluated with several tests. Tajima's D test yielded 
statistically insignificant results in all the four examinations 
(Korea, Tushima, Japan and all of East Asia) (р=0.084- 
0.755), whereas Fu's Fs test resulted in statistically signifi- 
cant results in all cases (р--0.00-0.009) (Table 4). Fu's Fs 
test showed high negative values (Korea, —24.78; Tushima, 
—6.41; Japan, —24.58; East Asia, —24.16), rejecting the 
assumption of neutral evolution of the control region locus 
or constant population size of P. altivelis over a long period 
of time (Table 4). 

The graphs resulting from the mismatch distribution analy- 
sis showed unimodal shapes in the groups, Korea and Japan, 
and the whole East Asian P. altivelis population data sets, 
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Table 3. AMOVA results from the hypothetical subdivision of seven groups (Korean mainland; JEJ, TUS, FKE; KYS, HNS, IWK, AMR, 
KZW, SAM, SAD; AMM; YAK; OKI; and China) based on hierarchical AMOVA (Appendix 1) 


Source of variation Degree of freedom 


Sum of squares 


Variance components Percentage of variation (96) 


Among groups 6 200.941 1.05312 Va 21.02 
Among populations within groups 14 82.185 0.17957 Vb 3.58 
Within populations 208 785.860 3.77870 Vc 75.40 
Total 228 1,068.987 


Fixation indices based оп 1000 permutations are like following: Фт=0.246 (p« 0.000), ®;=0.210 (p« 0.000) and Фс=0.045 (p< 0.006). 


The abbreviations used in Table 3 are listed in Table 1. 


Table 4. Neutrality test and mismatch distribution statistics using 244 mtDNA control region sequences 


Korea 


Tushima 


Japan East Asia 


Neutrality test 


Tajima's D —0.534 (p=0.340) 

Fu's Fs —24.782 (р-0.000)ғ 
Mismatch distribution 

SSD 0.00545 (pz0.490) 

HRI 0.00577 (p=0.890) 

t(tau) 10.488 (3.84-13.26)° 


SSD, sum of squared deviation; HRI, Harpending’s raggedness index. 
Significant at a=0.05, °99% percentile between lower and upper bound. 


respectively, except for Tushima group (Fig. 2). The null 
hypothesis of the sudden population expansion model could 
not be rejected in all groups either by the sum of squared 
deviation (p= 0.90) ог Harpending's raggedness index (p= 
0.85) (Table 4). 


DISCUSSION 


In the present study, we intended to use the phylogeography 
of the amphidromous ayu to understand the effects of paleo- 
ceanographic changes on marine organisms during the Qua- 
ternary period in East Asia. 

To get a more global picture of P. altivelis phylogeography, 
we substantially extended the sampling range of previous 
studies that have been conducted solely around the Japanesse 
archipelagos on the basis of mitochondrial DNA sequences 
(Iguchi et al., 1999; Iguchi and Nishida, 2000) to a large por- 
tion of East Asia, encompassing the Korean Peninsula. This 
resulted in a number of findings. First, population genetic 
structures and evidence of past demographic expansion of 
P. altivelis emerged more clearly as sampling range increases. 
However, when the analysis was confined to a specific local 
province, such as the Korean mainland populations, we could 
find no genetic structure, which suggests high gene flow rates 
among localities. Additionally, we could better understand a 
genetic connection among some islands in the Southern Sea 
of Korea, where island populations seem to be strongly influ- 
enced by the Tsushima Current. 
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0.532 (p=0.755) 
—6.411 (p=0.009) 


0.00406 (p=0.810) 
0.00681 (p=0.980) 
4.046 (0.56-15.12)» 


—1.293 (p=0.084) 
—24.582 (р=0.000)* 


—1.307 (p=0.062) 
—24.164 (p=0.002)? 


0.00041 (p=0.780) 
0.00477 (p=0.740) 
8.121 (6.18-10.89)* 


0.00052 (p=0.900) 
0.00312 (p=0.850) 
9.52 (5.81-13.01)* 


We should note that population genetic structures of P. 
altivelis became evident obvious as the examined area incre- 
ased. For instance, a significant differentiation in genetic 
structure was detected at the scope of East Asia (Фст--0.210) 
(Table 3), whereas the genetic pool of populations within 
the Korean mainland was shown to be highly homogeneous, 
with no significant genetic subdivision (Фст = —0.049). 
Moreover, we found that the genetic differentiation (high 
Фот values) in East Asian ayu has been caused mostly by 
distant insular, isolated lake and Chinese P. altivelis popula- 
tions (Table 2). Another factor that may have affected the 
patterns may be the amphidromous behavior of ayu, which 
migrate to the closest river rather than to their natal river 
during the marine stage (Kinoshita, 1993; Azuma et al., 2003; 
Iguchi et al., 2006). It has been suggested that the intensity 
of gene flow depends on the distance between populations 
by Iguchi and Nishida (2000). Additionally, a notably differ- 
entiated group was found to exist in three islands (JEJ, TUS, 
and FKE) located between the Korean peninsula and the 
Japanese archipelago. The most island populations tended to 
have high genetic differentiations to other populations, even 
among each island ones (Iguchi and Nishida, 2000). When 
these three islands were grouped together and then tested by 
the hierarchical AMOVA, it resulted in higher Фст values 
with the highest value of 0.210 (Appendix 1). Interestingly, 
these islands are located along the path of the Tsushima Cur- 
rent, which began to diverge from the Kuroshio Currents 
during the Quaternary period resulting in paleoceanographic 
changes in the region (Kitamura et al., 1999, 2001; Ujiié and 
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Ujiié, 1999). Thus, we can predict that the current might have 
influenced long-distance dispersal and genetic connections 
among the three insular P. altivelis populations. According to 
Teague et al. (2002), the dominant northeasterly Tsushima 
Current between JEJ and Tsushima Island (TUS) was esti- 
mated at about 25 cm/s across varying depths of sea water, 
from the surface to near the sea floor. If we apply this speed 
to the distance between JEJ and TUS (267 km), the duration 
for a passively neutral object from JEJ to reach TUS without 
considering friction or water mass is estimated at less than 13 
days. This duration may be enough for juvenile ayu to dis- 
perse to different islands within their marine life-stage before 
returning to rivers. Furthermore, the duration will shorten if 
the fish actively swim in the same direction as the current 
flows. Therefore, the high genetic homogeneity among the 
three insular populations strongly suggests a certain degree 
of realized migration between islands via the Tsushima Cur- 
rent during the ayu's marine stage. 

To understand the distribution and phylogeography of P. 
altivelis, it is important to consider the paleoceanographic 
changes in the region mentioned above. We can assume that 
cyclic sea-level and coastal line changes could have forced 
P. altivelis populations to experience a series of bottleneck 
and expansion events. To test this assumption, we examined 
the molecular DNA sequence data with neutrality test statis- 
tics, Tajima's D and Fu's Fs tests, and mismatch distribu- 
tion analysis. Fu's Fs statistics were significant in all data 
sets (Korea, Tushima, Japan and East Asia), suggesting either 
selection or sudden demographic expansion of the ayu popu- 
lation in recent past (Fu, 1997). The mtDNA control region, 
as a noncoding locus, is very unlikely to be under selection. 
Thus, the significantly negative values of Fu's Fs statistics 
likely indicate an expansion of population. The negative 
Tajima's D statistics may also indicate population expansion 
in the Japanese and East Asian data sets because the P-values 
of these statistics were close to 0.05 (0.084 and 0.062, respec- 
tively) (Table 4). Fu's Fs test can detect population expansion 
more powerfully than Tajima's D test (Fu, 1997). In addition, 
the unimodal shape and non-significant mismatch distribu- 
tion statistics (Fig. 2) support a sudden demographic expan- 
sion of P. altivelis in the recent past (Rogers, 1995). There- 
fore, we suggest that population demographic expansion and 
range expansion occurred simultaneously in ayu populations 
across East Asia. During the Pleistocene glaciations, temper- 
ature, sea level and costal line changes may also have great- 
ly affected both the demographic and distributional history 
of marine fishes (Lecomte et al., 2004; Atarhouch et al., 
2006). 

In parallel to the neutrality tests for the DNA sequences, 
the mismatch distribution analysis was used to address the 
question of recent population expansion in P. altivelis. Parti- 
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Frequency 


No. of pairwise differences 


Fig. 2. Mismatch distributions of mitochondrial control region 
sequences from Plecoglossus altivelis from different data sets; 
A, Korean mainland; B, Tushima group; C, Japanese mainland 
with Sado Island; and D, East Asian ayu group, respectively. 
Black line, observed value; Gray line, simulated value. 


cularly, the timing of population expansion (f, expansion 
time in generation) was estimated using the relationships: 
т--2иі and u=2uk (т, population expansion time scaled by 
mutation rate; и, mutation rate for the whole sequence; 4, 
mutation rate per nucleotide; and k, the length of the sequ- 
ence) (Excoffier et al., 2005). There is a range of divergence 
rates in the mtDNA control region across taxonomic groups 
of fish (3.6-2096) (Bowen and Grant, 1997; Donaldson and 
Wilson, 1999; Bowen et al., 2006; Ruzzante et al., 2008). 
The *time dependence of molecular rates" within shallow 
evolutionary time ranges (< 1 Myr) is important to consider 
(Ho et al., 2005; Penny, 2005; Ho and Larson, 2006). Thus, 
a range of rates of sequence divergence for galaxiid fish (G. 
vulgaris and G. divergens, 7.4-23.3% per million years) 
was chosen as a calibration point for molecular clock of the 
present study (Burridge et al., 2006). Assuming the genera- 
tion time as one year (Kim et al., 2005), the population expan- 
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sion for P. altivelis populations in East Asia (7=9.52) was 
estimated to begin since 126-391 thousand years (KY) ago. 
In the case of Korean P. altivelis, representing genetically 
homogeneous mainland populations, the expansion was esti- 
mated to begin about 138-430 KY ago. Likewise, P. alti- 
velis of the Japanese archipelagos began to expand around 
126-391 KY ago. These estimates correspond to the late 
Pleistocene and are sufficiently broad to encompass more 
than two glacial cycles before the LGM. Population expan- 
sion of other fishes during the late Pleistocene was also re- 
ported from the Yellow Sea to the East/Japan Sea. In the 
East China Sea and the Yellow Sea, the spotted tail goby, 
Synechogobius ommaturus (Song et al., 2010), that inhabited 
Asian inshore waters and one of demersal fish in Chinese 
Coastal waters, Nibea albiflor (Han et al., 2008), were esti- 
mated to begin expansion 84-248 and 85-170 KY ago, res- 
pectively. Also, the sailfin sandfish, Arctoscopus japonicas, 
that is distributed from the east coast of the Korea and the 
East/Japan Sea went through expansion about 90 KY ago 
(Shirai et al., 2006). In the case of fat greenling (Hexagram- 
mos otakii), the demographic history of population expan- 
sion was reported to occur during the late Pleistocene period 
(91-327 KY ago) (Habib et al., 2011). 

In addition, similar conclusions have been drawn with 
other marine pelagic fish. Liu et al. (2006) analyzed the his- 
torical demography of the Japanese anchovy, Engraulis japo- 
nicas, using mtDNA control region sequences and suggested 
79-317 KY ago as a population expansion time. Population 
of other pelagic fishes expanded during the last Pleistocene; 
Australian anchovy (E. australis, 45-178 KY ago) (Liu et al., 
2006), Pacific sardine (Sardina sagax, 229 KY ago), northern 
anchovy (Е. mordax, 311 KY ago) (Lecomte et al., 2004). 
Zhang et al. (2006) studied crimson snappers, Lutijanus ery- 
thropterus, and found genetic evidence of population expan- 
sion around 75 - 100 KY ago, which again supports the com- 
mon population dynamics of expansion around the late Pleis- 
tocene in East Asia. However, the range of expansion time 
in some species is too wide to match with a single glacial 
event in the past. Therefore it makes us difficult to pinpoint 
which past environmental factors really facilitated the demo- 
graphic increase in a wide array of marine fish species. The 
broad time estimates may be an artifact of the uncertainty in 
mutation rates of the DNA markers. We definitely need more 
accurate mutation rate estimates for the markers, although it 
is not likely to be easy and to come soon. Aside from the 
ambiguity of mutation rate in the control of mitochondrial 
DNA in marine fish, our estimate of expansion time since 
126-391 KY ago is surprisingly well before the LGM (19- 
26 KY) (Clark et al., 2009). This implies that even the severe 
climate change and sea level drop in the LGM could not 
effectively reduce the population size of P. altivelis to such 


Anim. Syst. Evol. Divers. 28(4), 279-290 


Populations Genetic Study of Ayu fish 


a degree that it might have erased the previous population 
genetic signature of growth. Again this interpretation intri- 
gues us what, then, really caused the population growth of 
ayu and other marine fish species since the middle to late 
Pleistocene period in the Northeast Asia. For a more compre- 
hensive phylogeography of P. altivelis in its full distribution- 
al range, more samples from China should be included, 
because the current data set does not represent the locality. 
Additionally it should be noted that multiple nuclear DNA 
markers are necessary for future studies to verify the present 
result based on only single mitochondrial locus data. 

In summary, the phylogeography of the ayu, P. altivelis, 
adds more evidence that there were prevalent population 
demographic changes in marine fish species characterized 
by population growth during the late Pleistocene long before 
the LGM in the Northeast Asia. Historical climate changes 
during the late Pleistocene, followed by sea level changes 
which would have inevitably led changes in the connectivity 
among marine populations, would have been common back- 
ground for the demographic change of ayu and many other 
marine fish. It is likely that the distribution and population 
genetic structures of these fish species vary depending on 
their life history characteristics and dispersal capabilities. 
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Appendix 1. Hierarchical AMOVA results based on mitochondria control region sequences from Plecoglossus altivelis in East Asia 
under various grouping methods 


Region groupings Фст value p of Фет Фс value p Of Dec 


(Korea, JEJ) (SAD, HNS, IWK, SAM, KZW, AMR, 


EE AHL AKC OKT ШЕКЕ (chine 0.09351 0.00098 0.15732 0.00000 
po Sie ee d ae ee ii 0.08494 0.00489 0.16465 0.00000 
,..... AMR; 0.05266 0.10557 0.18099 0.00000 
522. S pun PAPAE A 0.14348 0.00000 0.11482 0.00000 
2. c DUM E K 0.15566 0.00000 0.09973 0.00000 
polo SM KS. (ud Gene EE T 0.13425 0.00000 0.11179 0.00000 
VM NOW AMI s Chea 0.19328 0.00000 0.06888 0.00000 
TURON ME jn FERNER MAE 0.19734 0.00000 0.06517 0.00000 
prona nk ae an ey 0.20536 0.00000 0.05013 0.00000 
Шамал ATE CERIS ARMUT EMCO 0.21017 0.00000 0.04537 0.00684 
S IN CERA CMM CO, 0.14788 0.00000 0.09314 0.00000 
SK ARCANO CREAR AR RU. 0.19245 0.00000 0.04564 0.00000 
ped Noe US ERE TENEO AMR, 0.09351 0.00098 0.15732 0.00000 
pow pae qe aid DA M 0.08494 0.00489 0.16465 0.00000 
(Korea) (EJ, SAD; ANS; TWK; SAM: кай AMR, 0.05266 0.10557 0.18099 0.00000 


KYS, AMM, YAK, OKI, TUS, FKE) (China) 


Korean mainland population of P. altivelis labeled as “Korea”. 
Bold type, the highest Фет values. p, probability value obtained by permutation. 
The abbreviations used in Appendix 1 are listed in Table 1. 
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